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ABSTRACT 

Active Galactic Nuclei (AGN) play an important role in many aspects of the modern cosmol- 
ogy, and of particular interest is the issue of the interplay between AGN and their host galaxy. 
Using X-ray and optical data sets, we have explored the properties of a large sample of AGNs in 
the Subaru/XMM-Newton Deep Survey (SXDS) field, and studied their evolution in relation with 
the evolution of their host galaxy. We present here an anticorrelation between X-ray-to-optical 
flux ratio (X/O) and galaxy light concentration (C), which has been found for the first time and 
might suggest that early type galaxies, having poor matter supply to feed the AGN activity, have 
lower Eddington rates than those of late type galaxies. 

Subject headings: Galaxies: active — X-rays: galaxies 



1. Introduction 

Deep X-ray extragalactic surveys constitute an 
extremely useful tool to detect the population of 
active galactic nuclei (AGN), to investigate their 
time evolution, and to shed light on their trigger- 
ing mechanisms. Moreover, hard X-ray surveys 
are capable to detect all but the most absorbed 
(Compton thick) AGNs. Therefore, they provide 
the most complete and unbiased samples of AGNs. 
The combination of X-ray data with optical ob- 
servations allows succes fully tackling the stud 



of the AGN population (Georgakaki s et al 
Steffen et ahl 120061 : iBarcons et al.l 120071) 



•study 
2006; 



Sev- 



eral works dGebhardt et"al" 2000l ; Graham et al.. 
2001bl : iKauffmann et all 120031 )" have shown that 
the AGNs are directly related with some host 
galaxy properties, particularly with their bulges. 
A tight correlation between black hole (BH) mass 
and the b ulge dispersion velocity has been well es- 
tablished (jFerrarese fc MerritfeOOOHGebhardt et al 
2000), as well as with the galaxy light concentra- 
tion ([Graham et al.ll2001b[ ). 

We present here, for the first time, an anticorre- 
lation between the galaxy light concentration and 
the X-ray-to-optical flux ratio for a large sam- 
ple of AGNs in the Subaru/XMM Deep Survey 
(SXDS0 field, showing that more concentrated 
galaxies (earlier morphological types) have lower 
X/O ratio. This anticorrelation is found to be 
significant in a wide range of photometric red- 
shifts, with a similar first grade polynomial func- 
tion fitting in different redshift bins. We suggest 
that accretion rates in early type objects (having 
higher values of concentration index) in tipically 
gas-poor environment, are lower than in late type 
objects, that are surrounded by larger reserves of 
gas for AGN feeding. 

2. Observational data 

2.1. X ray data processing and source de- 
tection 

The XMM-Newton observations of SXDS field 



(jFurusawa et alJ |2008|) were collected from the 
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XMM-Newton Science Archive (XSA). Seven 
pointings around the central coordinates (a = 02^ 
18 m 00 s , <S = -05°00'00"; PI. M. Watson) have been 
gathered, covering a total area of ~ 1.3 deg 2 . 
The EPIC observations were reprocessed using 
SAS v7.1.2 standard reduction procedures. Six 
energy ranges were defined: soft (0.5-2 keV), hard- 
2 (2-4.5 keV), verydiard (4.5-10 keV), very.hard-2 
(4-7 keV), total (0.5-10 keV) and total-2 (0.5-7 
keV). Hardness ratios were defined as follows: 



HR (AiE, A 2 E) 



CR{A 1 E)~CR(A 2 E) 
CR{A 1 E) + CR(A 2 E) 



(1) 



where A\E y A 2 E are different energy bands and 
CR (A n E) is the count rate in a given energy 
band. The hardness ratio used in this work is 
HRh 2s =HR(hard-2, soft). 

Source detection was performed by means of the 
edetect_chain SAS procedure with a likehood 
parameter ML = — ln(l — P) > 14, where P 
is the probability that the source exists. We 
considered only those sources for which S/N > 
2 in the total band. We have detected 1121 
sources fulfilling these conditions. Limiting flux is 
1.3 x 10 -15 ergcm _2 s _1 in the 0.5-4.5 keV range. 

2.2. Optical counterparts 

Optical data have been gathered from the 
SXDS Public Suprime- Cam Data Release 1 (DR1) 



(jFurusawa et al.l [2008) . The total mapped area 
is 1.22 deg 2 , centered at a = 02H8"W, 5 = - 
05°00'00". Limiting AB magnitudes are 28.4, 
27.8, 27.7, 27.7 and 26.6 in B, V, R c , i' and z' 
filters, respectively. We have characterized the 
morphology of extended objects by means of the 
Abra ham concentration index (C; lAbraham et al 



1994J) , compu ted by a modified versio n of SEx- 
tractor 2.5.0 ( Bertin fc Arnoutsl Il996l modifica- 
tions from B. Holwerda), as the ratio between the 
integrated flux within a radius defined by the nor- 
malized radius a=0.3, and the total flux. 

We cross-matched the X-ray and optical cata- 
logs with a search radius of 3", obtaining a com- 
pleteness of 99.9% and a reliabilit y of 76.2% (ap- 



plying de Ruiter metho dology, see Ide Ruiter et al 



1977uPovic et al.l l2009. for more details). We have 
found 808 X-ray sources with optical counterpart. 
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3. Photometric redshifts and K-corrections 



sented in Fig. Q] (already observed in the popula- 



We have used ZEBRA (jFeldmann et al.l l2006h 



tion of AGNs of the Groth field, see Povic et al 



to compute photometric redshifts of the optical 
counterparts. Input photometric data includes 11 
optical and IR bands: Subaru/SuprimeCam B, 
V, R c , i' and z'; J, H and K from the UKIDSS^ 
survey; Spitzcr/IRA C DR2 data in 3.6^tm 4.5/im 
and 5.8/im bands ( Suracd 120051 ). The best 
SED template set has bee n gathered from th e 



SWIRE template library (jPolletta et all l2007f) . 



As a complementary check, photometric redshifts 
have been also comput ed using the HyperZ code 
(|Bolzonella et al.l [2000) with the same template 
set, accepting only those objects for which the 
redshift differences from both codes are less than 
0.1 for z<l and 0.2 for z>l. From the initial 
sample of 808 sources with optical counterparts, 
377 objects fulfill these conditions. 

In order to obtain K-corrections of the opti- 
cal fluxes we have used the idl routine KCOR- 

On the other hand, 



2007) 



rect (jBlanton et al. 
the X-ray fluxes were K-corrected by assum- 
ing a standard power law SED with T = 1.8. A 
rough classification between unabsorbed and ab- 
sorbed AGNs has been made u sing HRft 2si that 



is qu it e sensitive to absor ption (jDella Ceca et al 



2004t iDwellv et all l2005h : those objects with 
HRf l2s < 0.35 (some 53%) have been considered 
as unabsorbed, and no intrinsic absorption has 
been applied to the power law SED; otherwise 
the objects have been considered as absorbed and 
a fixed intrinsic absorption Nh = 1-0 x 10 22 cm~ 2 
has been included. 

4. The anticorrelation between galaxy 
light concentration and X/O flux ratio 

Correlating the X-ray and optical properties we 
found a very strong anticorrelation between the 
X/O ratici and the concentration index, repre- 



2 The UKIDSS project is defined in iLawrencd J2007f> 
UKIDSS uses the UKIRT Wide Field Camer a (WFCAM) 



The photometric system is described in Hcwett ct al 
(2006), and the calibration is described in iHodgk in et al 
(2009). The pipeline processing and scien ce archive are 
descri bed in Irwin et al (2009, in prep) and Hambly ct al. 
(2008). We have used d ata from the 3th data release, which 
is described in detail in I Warren et al.l (120081 ). 

3 computed as i^o.S— 4.5kev/^ ? Hi where the optical flux Fr 
has been derived from the SExtractor "auto" magnitudes 
in the R band 



2009t ISanchez-Portal et all 1200a ). The X/O ra- 
tio measures the X-ray flux (in 0.5-4.5keV band), 
normalized with the R-band optical flux of the 
whole galaxy (nucleus, bulge and disc), although 
the sources in our sample are predominantly nu- 
cleus/bulge dominated. For the whole sample, 
a first order polynomial function has been fitted 
with a slope of -2.76 ±0.14 and log X/O-intercept 
of -1.54 ±0.07. Both unabsorbed and absorbed 
AGNs follow the same relation. Final C and X/O 
error bars are shown in the first plot of Fig. [TJ and 
are significant in all six redshift ranges. In general, 
we obtained that C is accurate to 5%, mainly due 
to the uncertainty introduced by the nuclear-point 
source influence (see Section |4~T|) . The average er- 
ror of X/O flux ratio is AX/O = 0.351 for the full 
sample. 

4.1. Analysis of possible bias effects 

- Redshift and Bandpass shift. It is expected 
that distance will directly affect the measure- 
ments of C, tending to shift objects to appar- 
ently later Hubble ty pes ( Brinchmann et al. 19981 : 
Conselice et al.ll2003l ): besides this, there is an in- 
fluence of a bandpass shift with res pect to the rest 
wavelength ( Bershady et al. 2000l ). We obtained 
K-corrected X/O ratios and we created six red- 
shift ranges, as shown in Fig. [TJ The range of C 
(and X/O) is the same in all six bins, i.e. basically 
independent of redshift. Moreover, the relation is 
the same in all redshift ranges, as verified by linear 
regression fits. These facts indicate that the effect 
of these biases is negligible. 

- Obscuration. Intrinsic absorption of the nuclear 
point-source could make the X/O flux ratio larger 
and C smaller, leading to apparent anticorrelation 
between these two parameters. We used HR/ l2s 
criterion in order to examine the absorption effect 
on our anticorrelation. As already mentioned in 
Section [3l this parameter should select efficiently 
the majority of the most heavily absorbed sources. 
Fig. [2] shows the relationship between HR^s and 
C in the same six redshift ranges as in Fig. [TJ No 
correlation was found in any redshift range. So, 
even if HRj,2 S might not be effective to discrim- 
inate absorbed AGNs at high redshifts, we can 
at least claim that for objects within z< 1.5 (e.g. 
Akylas et al.ll2006l ) obscuration should not play an 
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Fig. 1. — Relationship between the concentration index and the K-corrected X/O ratio in six redshift bins. The 
continuum line corresponds to the linear fit for data points in each redshift bin, while the dashed line represents the 
fit for the whole sample. Final error bars are included in the first plot, and are significant in all six redshift ranges. 
We found that C is accurate to ~5%. The average error of X/O ratio is found to be AX/O = 0.351 for the full sample. 
Goodness of full sample fit, as measured by the Spearman's rank correlation coefficient (p) is -0.708. 



important role in the observed anticorrelation. 
- Nuclear point source. The presence of a nuclear 
point-like source could affect the determination of 
C, in the sense that the measured values are over- 
estimated for the host gala x y com paring with their 
true values. Gabor et al. ( 2009t ) found that the 
presence of a nuclear point source could have a 
significant bias on concentration measurements to 
higher values, although the authors do not quan- 
tify it. However, the influence of the nuclear point 
source in the concentration indexes has been inves- 
ti gated in HST count erparts of Chandra sources 
bv lGrogin et al.l (|2003h . finding that the AGN host 
concentration differences are truly related to dif- 
ferences in the host galaxy structure and not to the 
nuclear point-source flux contribution, quantify- 
ing the possible influence at C in ~ 10%. We have 
extended this analysis to ground-based data, by 
building several galaxy models (elliptical, lentic- 



ular, early, intermediate and late spirals), con- 
volved with an HST-like PSF (a gaussian func- 
tion with FWHM = 0.043 arcsec) and a ground- 
like PSF (FWHM = 1.0 arcsec) representative of 
our data, adding in both cases different fractions of 
an unresolved nuclear source (nucleus/total frac- 
tion up to 50%) and deriving the corresponding 
Abraham concentration indexes. Even though the 
range of C of ground- like data shrinks considerably 
with respect to that observed in HST-like data, 
we observe that HST and ground concentration 
indexes are very tightly linearly correlated and, 
that the variation of C is dominated by the host 
galaxy morphology and not by the nuclear compo- 
nent contribution. Both results suggest that the 
influence of this bias in the anticorrelation is small. 
We derive that the maximum effect of the nuclear 
point source on C is about 7% in HST-like models 
and less than 3% in ground-like ones. 
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- Seeing, Apparent Brightness and Size. Optically 
fainter objects tend to show lower concentration 
indexes, and this effect is enhanced with poorer 
seeing. In fact, this correlation is observed in our 
data: the apparent magnitude and logC are lin- 
early correlated^- A weak tendency that appar- 
ently smaller objects have lower concentration in- 
dexes is present in our data, affecting mainly late 
type galaxies. Besides the observational bias ef- 
fect, this tendency could be also partially intrinsic 
(as in the case of the apparent brightness), show- 
ing that intrinsically smaller objects have lower 
C. In order to test if C is still a good tracer of the 
host galaxy morphology, we have performed an in- 
dep endent classificati on based on visual inspection 
(see lPovic et al. 2009). There is a good agreement 
between the results derived from the combination 
of the concentration and asymmetry indexes and 
those obtained by inspection. Therefore, even tak- 
ing into account that apparent brightness and see- 
ing effects play an important role, we are still con- 
fident that C is a good tracer of galaxy morphol- 
ogy, and that, on average, lower concentration in- 
dexes are tracing later morphological types. More- 
over, we performed a preliminary ana lysis of pub- 



lic HST/ACS images o f the EGS field (|Davis et al 



20071 ; lLotz et al.l [2008) , obtaining the I-band op- 
tical fluxes and concentrations; we combined the 
optical infor mation with the A EGIS EGS X-ray 
public data ( Laird et al. 20091 ). obtaining a clear 
anticorrelation between the X/O ratio and C, in 
all four AEGIS X-ray energy bands. 

4.2. Comparison with galaxy mock 
catalogs 

We used The Millennium- II Simulations database 
( Bovlan-Kolchin et al.|[2009l . and references therein) 
in order to perform further tests of the existence 
of an anticorrelation between X/O and C. The 
mock catalog we used, contained in the MPA- 
Galaxies Database, is the table DeLucia2006a 
(jCroton et all 120061 : 



de Lucia fc Blaizotl l2007h 



We worked with more than 700000 galaxies, divid- 
ing them in the same redshift ranges as in Fig. [T] 
We obtained the R-band bulge-to-total flux ratio 
(B/T) for all simulated galaxies; this parameter 



4 Actually this correlation is twofold: besides the observa- 
tional bias discussed here, it is observed that less con- 
centrated galaxies tend to be intrinsically fainter (e.g. 
iBlanton et al.ll200ir ): this trend is also found in our data. 



is proportional to the galaxy light concentration. 
The X/O ratios are calculated from the catalog's 
X-ray bolometric luminosities and the R-band 
bulge fluxes. Fig. [3] represents the obtained rela- 
tionship for more than 60000 simulated galaxies 
with z < 0.4. A clear anticorrelation can be seen 
for all B/T values except for pure bulge galaxies 
(i.e. with log(B/T) > -0.25), whose real popula- 
tion should be very small, so we do not expect 
to have a significant number of these galaxies in 
our real sample. In the log(B/T) interval « [- 
1.5,-0.2], where most of our objects are located, 
this clear anticorrelation can be confirmed, show- 
ing that less-concentrated or less bulge-dominated 
objects (later Hubble types) have higher X/Obuige 
flux ratios. Finally, for galaxies with log(B/T) <- 
1.5, we can expect to be below our X-ray de- 
tection limits. There is a significant number of 
outliers, mostly found below the observed anti- 
correlation. All these objects show properties of 
small, dwarf or satellite galaxies (small maximum 
rotational velocitiy, reduced virial and stellar mass 
and hot/cold gas content, low SFR) and low bolo- 
metric X-ray luminosities, well below 10 42 ergs" 1 , 
making most of them undetectable at our sensitiv- 
ity thresholds. In our sample, only 3% of all X-ray 
emitters with optical counterparts have X-ray lu- 
minosities lower than I0 42 ergs _1 , in the observed 
0.5-4.5 keV energy range. Figure [3] shows a clear 
difference in X-ray luminosity between the galax- 
ies belonging to the anticorrelation sequence, and 
to the outliers. However, there are examples of 
outliers found observationally (see Section [5]). 

5. Discussion 

All the tests performed, described in Sections 
14.11 and 14.21 suggest that the observed anticorre- 
lation cannot be solely attributed to biases, but 
also to a true physical effect. Deriving a possible 
explanation requires understanding the physical 
connections between the observables C and X/O. 
There is a tight correlation between the C31 
concentration index (reflecting morphologies of 
bulges) and the mass of the central supermas- 
sive BH, showing that m ore concentrated galax ies 
have more massive BHs ([Graham et al.lEoOf btf al). 
The Abraham concentration index C used in this 
work is not the most effective tracer of bulge mor- 
phologies; however, all the tests performed so far 
suggest that C is a parameter good enough to 
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Fig. 2. — Relationship between the concentration index C and the HR^s hardness ratio in the same redshift bins 
as in Fig. [T] 



trace the galaxy light concentration and that, on 
average, higher concentration indexes are tracing 
more concentrated objects (or earlier morpholog- 
ical types), while lower concentration indexes are 
tracing later morphological types. Moreover, C 
and C31 concentration index are monotonically 
correlated. We analyzed the morphology of SXDS 
X-ray emit ters with optical counterparts using 
GALSVM (iHuertas-Companv et all l2008lh We 
obtained Gini and Moment of light (M20), finding 
a clear correlation between these parameters and 
C, indicating again that higher values of C are 
tracing earlier Hubble types. On the other hand, 
there is a tight correlation between Gi ni and M2n 



(.e.g. 



Lotz et al 



with C31 concentration index 
2004]) . 

On the other hand, the physical connections of the 
X/O flux ratio are not so evident, but as shown 
below it could be interpreted as an indicator of 
the Eddington ratio. As already mentioned, our 
sample is predominantly made of nucleus/bulge 
dominated sources. Thus we can distinguish be- 



tween two possible cases: 

- If the nuclear luminosity dominates that of the 
host galaxy, the X/O ratio can be thought as a 
measure of the X-ray to optical spectral index, 
a o jJE On the other hand, there is a strong corre- 
lation between the hard/soft X-ray s pectral index , 
cex, and the Eddington ratio (e.g. lGrupej|2004 ; 
lBian| [2005h . The relationship between ax and 
aox has not been thoroughly studied in large 
samples of AGN, but a linear correlation has been 
found between both spectral indexes in a sample 
comprising 22 out of 23 quasars in the complete 
the Palom ar-Green X-ray sam ple with z < 0.4 and 
M B <-23 (|Shang et al.ll2007l ). If this last correla- 
tion holds for our sample, X/O could be tracing 
the Eddington ratio in the large nuclear luminos- 
ity limit. 

- If the host galaxy luminosity is large when com- 
pared with the nuclear luminosity (Lbuige+disc 3> 



5 This relation is usually defined as ctpx 
0.3838 log(/ 2feeV // 25ooA ) 
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Fig. 3. — Relationship between t he X-ray-to-bulge optical flux ratio and bulge-to-total flux ratio (R band) for The 
Millennium-II simulated galaxies (|de Lucia fc Bl aizot 2007) in the redshift range z < 0.4. The same relationship is 
obtained in the other five redshift ranges shown in Fig. [T] 99% of the objects located below the dashed line have 
X-ray bolometric luminosities less than 10 ergs" 1 , i.e. below the faint end of our sample. 



Lnudeus), the X/O ratio can be thought as a lower 
limit of the X-ray-to-bulge luminosity ratio, that 
is in turn a (weak) measure of the AGN Eddington 
ratio L/LEdd °c L/Mbh, assuming that the X-ray 
luminosity represents the nuclear luminosity and 
the bulge luminosity is proportional to the bulge 
mass (and therefore to M^). 
Hence, we can guess a correlation between X/O 
ratio and the energy production efficiency of the 
AGN measured by the Eddington ratio. Under 
these assumptions, the galaxy light concentra- 
tion vs. X/O flux ratio relation traces the cor- 
relation between the nuclear BH mass and the 
Eddington ratio. This result could therefore sug- 
gest that more concentrated or bulge-dominated 
(early- type) galaxies, having poor matter supply 
to feed the activity, have lower Eddington rates 
than those of late-type, with larger reservoires of 
the material for AGN feeding. 



This suggested ap proach is con s istent with the 
results obtained by iBallo et alJ (|2007l ). They 
found that AGNs with large SMBH (M B h > 3 x 10 6 M©), 
have low X-ray luminosity and Eddington rate <C 
1 , and conversely, that smaller SMBH have higher 
luminosity in X-rays and Eddington rate ~ 1, 
which corresponds to a more efficient accretion 
rate. 

O ur approach is al so consistent with the results 
of lWu fc Liu] (|2004l ). which studied the Mbh and 
Eddington rates of a sample of 135 double-peaked 
broad line AGNs. They obtained that if the sep- 
aration between the line peaks decrease (this sep- 
aration is directly related with the FWHM of the 
line, and correlated with the Mbh), the Edding- 
ton rate increase s. 

Kawakatu et al.l (|2007| ) found an anticorrelation 



between the mass of a SMBH and the infrared-to- 
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Eddington luminosity ratio, L/^/L^d, for a sam- 
ple of type 1 ULIRGs and nearby QSOs, which 
again could be consistent with our approach. The 
anticorrelation is interpreted as a link between the 
mass of a SMBH and the rate of mass accretion 
onto a SMBH, normalized by the AGN Edding- 
ton rate, which indicates that the growth of the 
BH is determinated by the external mass sup- 
ply process, and not the AGN Eddington-limited 
mechanism, changing its mass accretion rate from 
super-Eddington to sub-Eddington. 
Nevertheless, as pointed out in Section 14.21 there 
are observational examples of outliers of this anti- 
correlation, showing that massive elliptical galax- 
ies can have high accreting BHs, and that late- 
type galaxies can host weakly accreting engine s 
(e.g. Churazov et al. 20051 : Desroches et al. 2009). 
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